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ABSTRACT 


Previous  analytical  studies  endeavoring  to  accomt  for 
vortex  induced  rolling  nxDments  on  missiles  have  not  successfully 
correlated  with  e:q5erimental  results.  This  points  up  the  need  for 
a  theoretically  more  sophisticated  treatment,  which  the  present 
work  attempts  to  satisfy. 

A  model  for  the  two  dimensional  cross  flow  past  the  tail 
body  combination  is  employed  and  the  Blasius  complex  contour 
integral  used  to  find  the  moments.  The  integral  is  solved  by 
examining  the  singularities  (poles  and  branch  points)  inside  the 
contour. 

This  is  the  first  scientific  report  in  a  continuing 
research  effort.  The  final  result  contained  herein  presents  the 
complete  mathematical  solution  for  the  two  dimensional  ideal  flow 
past  a  planar  tail  body  combination.  It  is  in  closed  analytic 
form  with  the  exception  of  one  definite  integral  which  lias  been 
reduced  to  a  form  suitable  for  numerical  calculation.  All  terms 
of  the  solution  (integrated  as  well  as  the  integral)  are  best 
calculated  for  various  values  of  the  parameters  on  a  large  scale 
digital  computer. 
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SECTICTJ  !♦  Introduction  and  Review  of  Previous  Work 

The  successful  research  and  development  leading  to  modem 
rocket  flight  has  brought  with  it  several  secondary  aerodynamic 
problems  to  be  solved  as  a  result  of  the  unstable  behavior  that 
these  rockets  occasionally  exliibit.  Some  types  of  instability 
that  have  been  ejgjerienced,  in  particular  by  sounding  rockets, 
is  caused  by  anomalous  rolling  moments  which  thus  far  have  not 
been  successfully  predicted  by  theory. 

Physical  misalignments  of  the  vehicle  will  tend  to  cai:ise 
lack  of  symmetry  about  the  roll  axis.  Hence  the  vehicle  is  usually 
made  to  spin  by  canting  the  fins,  in  order  to  minimi25e  any  non- 
syirenetric  behavior  during  flight  through  the  atmosphere.  In 
addition,  spin  provides  a  gyroscopic  type  stability  in  the  endo- 
atmospheric  environment.  The  roll  rate  is  usually  proportional  to 
the  rocket  free  stream  velocity.  However,  as  described  in  reference 
1,  disastrous  results  often  occur  when  the  roll  rates  deviate 
significantly  from  the  expected  values.  For  exanple,  if  the  roll 
rate  becomes  the  same  as  the  natioral  pitching  frequency  of  the 
vehicle,  thei.  pitch-roll  coupling  is  liable  to  occur.  That  is,  the 
rocket  experiences  the  maximum  angle  of  attack  in  the  same  plane; 
any  misalignments  tend  to  increase  the  angle  of  attack  causing  high 
drag  in  the  atmosphere  and  high  nutational  movements  outside  of  it. 

Another  type  of  instability  is  caused  by  the  Magnus  effect; 
due  to  spin  there  is  a  lateral  force  on  the  cylindrical  body  of  the 
vehicle  in  a  direction  which  is  perpendicular  to  the  angle  of  attack  plane 
(i.e.  the  plane  consisting  of  the  free  stream  velocity  vector  and 
the  body  roll  axis).  (Ref.  2.)  It  may  be  shown  that  the  roll  induced 
force  on  the  fins  is  in  the  opposite  direction.  (Ref,  3.)  Thus, 
since  the  centers  of  pressure  for  the  body  and  fins  are  separated, 
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there  is  a  yawing  moment  on  the  vehicle.  Needless  to  say,  unpre¬ 
dictable  roll  rates  (and  therefore,  to  a  large  extent,  mcontxol- 
lable  roll  rates)  tend  to  aggrevate  flight  stability. 

Many  strange  types  of  spin  behavior  have  been  experienced  by 
reentry  vehicles  and  torpedoes  as  well  as  rockets.  These  include 
spin  reversal  and  constant  velocity  spin,  even  in  cases  whiere 
canted  fins  are  part  of  the  configuration.  There  have  been  several 
attempts  in  the  literature  to  explain  the  unexpected  rolling  moments. 
For  example,  thermo-elastic  effects  on  the  fins,  (Ref.  4),  and 
oblique  heat  shield  tape  winding  on  the  noses  of  reentry  vehicles 
(Ref.  5)  have  accounted  for  small  scale  behavior,  but  have  been 
heretofore  unsuccessful  in  explaining  some  of  the  large  moments 
that  are  obseirved.  In  particular,  there  is  need  for  analytical 
wc~k  in  the  low  roll  case,  as  for  sounding  rockets  v/here 

is  the  l-ingential  velocity  of  the  cylindrical  surface  due  to 
roll ,  and'\^  is  the  free  •  _ream  velocity. ) 

One  approach  that  shows  promise  of  explaining  the  rolling 
moment  phenomena  in  the  low  roll  case  is  as  follows:  ai  high 
angle  of  attack  vortices  are  formed  at  the  rear  of  the  cylindrical 
body.  When  these  vortices  are  shed,  the  longitudinal  free  stream 
velocity  (the  component  along  the  body  axis)  washes  them  past  the 
fins  resulting  in  inteference  that  affects  the  rolling  moment. 

Figure  1,  which  has  been  taken  from  Ref.  6,  depicts  the  manner 
in  which  the  vortices  are  foirraed  on  the  non-finned  part  of  the 
vehicle.  Figure  2  (a),  taken  from  Ref.  7,  is  a  photograph  in  the 
cross  flow  plane.  These  vortices  are  symmetrical;  however  a  roll 
velocity  will  cause  non-syinmetrical  vortices  as  shown  in  Fig.  2  (b) 
(Ref.  7)  for  which  17^  j  . 


Several  attempts  have  been  made  in  the  literature  to  deal 
with  vortex  induced  moments,  but  the  analytical  methods  used  have 
not  been  exact,  and  largely  unsuccessful  in  predicting  e:q)erimental 
results. 


One  approach  (Ref.  6)  uses  the  strip  theory  method:  firstly 
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the  flow  in  the  cross  plane  is  considered  inviscid  and  incom¬ 
pressible  (see  Section  II  for  a  justification).  Then  potential 
fliiid  theory  is  employed  to  find  the  flow  past  a  circular  cylinder 
only,  due  to  an  external  vortex,  P  ,  by  suitably  placing  an 
image  vortex  inside  the  cylindcT.  A  local  angle  of  attack,  , 
on  a  fin  as  a  function  of  the  spanwise  distance  y  (see  Pig.  3) 
is  described  by  dividrng  the  velocity  induced  by  the  vortex  normal 


Fig.  3 

to  the  fin  surface,  1/^  ,  by  the  free  stream  velocity,  ^ 
Then  the  element  of  normal  force  on  the  fin  is  taken  to  be 

cIM  = 

where  C  is  the  chord  length,  Q  the  dynamic  pressure  and, 
assximing  linear  theory,  ^ 

r  =  ^ 

c<. 


o 


is  the  normal  force  coefficient  slope.  The  rolling  moment  on  one 
fin  is  given  by 


X  ci 


and  thus 


J5!_  ;<  C  (f.)  ciy  , 


The  induced  rolling  moment  on  a  multi-fin  configuration  is  then 
obtained  by  addition.* 

Rolling  moments  calculated  by  the  preceeding  method  using 
symmetric  vortices  did  not  agree  with  those  obtained  from  wind 
tunnel  studies.  (Ref.  6,  p.  35)  Obviously  strip  theory,  although 
successful  in  some  applications,  has  many  shortcomings.  First  of 
all  the  streamlines  due  to  a  vortex  in  the  presence  of  a  cylinder 
only,  is  not  the  same  if  a  fin  is  added;  indeed,  by  the  time  the 
vortex  reaches  the  fin,  the  streamlines  are  distorted  so  that  there 


is  no  normal  velocity  component.  To  suppose  that  the  normal  force 
curve  slope, ,  on  a  flat  p3,ate  caused  by  a  tr^ve  free  stream 
angle  of  attack  is  the  same  as  that  due  to  theo^  defined  above, 
is  an  approximation  at  best.  Certainly  a  more  accurate  model  would 
be  to  find  the  vortex  flow  in  the  presence  of  the  fins  as  well  as 
the  cylinder.  Secondly,  even  if  the  moment  on  one  fin  could  be 
correctly  computed,  the  total  moment  due  to  all  the  fins  cannot 
be  correctly  obtained  by  addition,  since  each  fin  interferes  with 


•The  appropriate  flat  plate  coefficient  for  supersonic  linear 
theory  is  (Ref.  8)  ^  14- 


—  — 


7 


and  obstructs  the  flow  past  the  others.  While  strip  theory  might 
give  good  predictions  of  moment  for  a  single  fin  in  the  vicinity 
of  the  vortex,  multiple  fins  present  a  much  too  complicated  flow 
field  to  be  attacked  by  the  simple  model.  Thirdly,  the  effect 
of  the  con^jonent  of  free  stream  velocity  in  the  cross  flow  plane 
is  not  accoiaited  for,  nor  is  a  possible  Kutta  condition  in  this 
plane  considered. 

Reference  9  is  a  continuation  of  the  vrork  of  Ref.  6, 
and  attempts  to  improve  results.  It  uses  a  better  expression 
for  the  normal  force  coefficient  by  taking  account  of  non¬ 
linear  hypersonic  terms;  the  expression  is  equivalent  to 

r  ^  d£jL  ^ -i- dl£^  -i- .  ^  ^ 

However,  comparisons  between  e^qjerimental  data  and  theory  showed 
no  improvement.  (Ref.  9,  p.  18.) 

In  Ref.  10  still  another  theoretical  approach  is  presented 
slender  body  theory  is  used  in  conjunction  with  the  momentum  method 
to  get  an  order  of  magnitude  estimate  for  rolling  moments.  However 
one  needs  experimental  values  of  the  base  pressures  in  order  to 
apply  the  momentum  method. 

There  have  been  attempts  to  find  vertex  induced  moments 
experienced  bv  varying  cross  sectional  slender  airplanes  in 
unsteady  flow  using  the  Blasius  method.  (Ref.  11  and  12. )  However, 
the  correctness  of  this  work  is  questionable.  One  of  the  short¬ 
comings,  among  others,  is  the  failure  to  account  for  the  fact  that 
Blasius'  formula  does  not  allow  for  varying  cross-section  in 
unsteady  flow.  (See  Ref  13,  p.  119.) 


The  general  model  which  we  would  like  to  consider  con¬ 
sists  of  inviscid  flew  past  a  long  slender  cylindrical  body  as 
depicted  in  Fig.  1,  with  a  tail  configviration  as  shown  in  Fig.  4. 
The  cant  of  the  fins  is  considered  negligible,  and  the  roll  rate 
of  the  vehicle  so  low  that  steady  state  flow  may  be  assumed.* 


Now  consider  the  flow  past  the  vehicle  shown  in  Fig.  5.* 
It  is  assumed  initially  that  the  fins  are  rectangvilar  flat  plates 


*  The  rolling  moment  to  be  found  usi;.g  these  assumptions  will  be 
added  to  that  due  to  fin  cant,  assuming  that  linearity  allows 
for  superposition. 

‘•The  effective  angle  of  av.tack  is  taken  to  include  sideslip 
velocity. 


(although  the  assumption  of  unvarying  span  can  be  relaxed  as  will 
be  described  in  Section  IV).  The  free  stream  velocity  may  be 
divided  into  the  longitudinal  component, 

cross  flow  velocity  Vao  ^  .  Because  the  fins  are  symmetrically 
placed  and  negligible  roll  is  assumed,  the  longitudinal  velocity 
component  should  have  little  effect  on  the  rolling  moment  for  non- 
canted  fins.  The  force  on  any  fin  due  to  this  component  is 


Figure  5 

negated  by  the  force  on  the  opposing  fin,  180°  away.  (Since 
the  longitudinal  flow  is  basically  paraxlel  to  the  long  cylindrical 
surface  it  is  unobstructed  until  it  reaches  each  fin,  and  so  should 
have  the  same  effect  on  all.)* 

*  Naturally,  the  longitudinal  flow  will  induce  a  moment  on  differen¬ 
tially  canted  fins. 


10 


Or  the  other  hand  the  cross  flow  velocity  is  a  majcr  factor  in 
rolling  moments.  First  of  all  it  is  the  cause  of  the  vortex 
generation  and  shedding  depicted  in  Pig.  It.  Secondly  the  cross 
flow  shou3.d  affect  opposing  fins  differently  since  the  cylindrical 
body  presents  an  obstructive  infliience  as  the  fliaid  goes  past  it. 

The  main  role  of  the  longitudinal  velocity  is  to  wash  the  vortices 
produced  on  the  ma?-n  body  siu:face  aft,  past  the  tail  body  ccmbin- 
ation. 

For  a  typical  free  stretm  Mach  number  of  5  and  an  angle 
of  attack  egival  to  20*,  the  cross  flow  Mach  n'unu;er  is  about  1.7. 
However,  when  we  take  the  conical  bow  shock  into  consideration, 
(whose  geometry  is  determined  by  the  vehicle  nose  configuration), 
the  cross  flow  velocity  is  further  reduced,  and  must  be  subsonic 
in  at  least  a  considerable  region  behind  the  shock. 

If  one  were  interested  in  yawing  or  pitching  !-ioments, 
the  entire  body  of  the  vehicle  would  play  a  role  in  the  force 
distribution.  Hence,  since  the  vehicle  contour  changes  in  the 
longitudinal  direction,  a  three  dimensional  treatment  (usually 
approximated  by  slender  body  theory)  would  be  called  for.  In  the 
case  of  roiling  moment,  hov/ever,  it  is  only  the  fins  which  e;:pe''- 
ience  contributing  forces  in  inviscid  flow,  and  since  they  are 
rectangular,  the  two  diriensional  cross  flow  problem  will  serve  as 
oxxc  model.  Because  this  flow  has  been  shown  to  be  largely  s\±>- 
sonic,  the  approximation  will  be  made  that  the  cross  flow  Mach 
number  is  zero.  It  may  later  be  necessary  to  use  perturbation 
methods  to  take  account  of  compressibility  effects. 

The  final  decision  in  describing  the  model  is  to  choose 
the  tail  body  configuration.  Although  the  crucifcrm  shape  is  common 
to  the  sounding  rocket  prototype,  the  ideal  fluid  which  we  have 
assumed  would  be  questionable  because  flow  separation  is  likely  to 
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occur  for  any  roll  angle  (p  (see  Pig.  4  (a)).* **  Moreover,  the 
separated  flow  region  is  likely  to  be  in  the  important  vicinity 
of  the  shed  vortices.  Since  a  satisfactory  mathematical  model 
of  separated  flow  phenomenon  is  lacking,  we  will  choose  our  first 
ph-ioical  model  to  be  the  planar  tail  body  configuration  (Fig,  4(b)). 
Then,  at  least  for  small  roll  angles,  the  question  of  separation  is 
not  as  serious.  Moreover,  this  model  lends  itself  to  taking  account 
of  compressibility  effects  by  making  a  Prandtl-Glauert  transforma¬ 
tion,  since  the  planar  cross  section  obviously  better  satisfies 
the  slender  body  requirement  (for  small  )  than  does  the  cruci¬ 
form  configiuration.  Although  the  planar  tail  body  does  not  rep¬ 
resent  the  usual  prototype  it  is  a  simplier  mathematical  model  and 
one  which  is  more  consistent  with  the  two  dimensional  physical 
reality.  Our  main  purpose  is  to  obtain  an  understadning  of  the 
fundamental  phenomena  underlying  vortex  induced  moments,  (fjce  this 
understanding  has  been  accomplished  and  successful  correlations  made 
with  experiment,  then  it  would  be  in  order  to  attempt  to  formulate 
a  model  for  flow  past  a  cruciform  tail  as  the  next  step. 

Let  ^  be  the  cross  flow  velocity,  s<  ,  past  the 

configuration  shown  in  fig.  4  (b),  reduced  by  an  appropriate  factor 
due  to  the  conical  bow  shock  v/ave.  ••  The  roll  angle  ^  is  the  degree- 
to  which  the  fins  have  turned  out  of  the  vehicle's  angle  of  attack 
plane-.  Let  us  represent  the  complex  potential  for  this  flow  in  the 
2  plane  by  If  (2)  (see  Fiy,  6}  where  acco’cit  will  be  taken  of  two. 


*  The  concern  over  separation  i<;  in  a  laboratory  two  dimensional 
experiment.  In  the  rocket  prototype,  the  longitudiMO..'  velocity 
component  would  largely  eliminate  any  cross  flow  separation  effects. 

**T7  is  c'ssumc-d  to  be  constant  since  the  bow  shock  can  be  taken  to 
be  revHsonably  flat  in  the  vicinity  of  the  tail  section. 


Figure  6.  The  ^  plane 


not  necessarily  symmetrical  vortices ,  Pj  and  /^  •  *  Then  the 
steady  state  rolling  mor>ent  about  the  origin  is  given  by  the 
Blasius  aquation 

<i)  2: 

where  the  contour  is  along  the  outside  of  the  body  surface. 

(Ref.  13,  p.  49.) 

It  will  be  most  e^qjedient  to  evaluate  Equ.  (1)  by 
mapping  the  body  contour  into  a  circle.  This  can  be  accomplished 
by  a  two  step  Jovikowsky  transformation: 

-a:'  - 

(2)  £  - 

which  maps  the  cross  section  into  a  flat  plate  (Fig.  7) 

•  It  is  clear  that  the  strip  theory  method  for  finding  the  moments 
is  no  longer  appropriate  here  since  the  velocity  represented  by 

=:  ciir/d  i  has  zero  normal  component  to  the  fin  surface. 
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and 


(3) 


s 


•+- 


c 


2- 


which  maps  the  flat  plate  into  a  circle,  (Fig.  8). 


Figure  7,  The  g'  plane 


The  points  c  and  -c  in  the  ^  plans  are  mappings  of  '6~ 
and  in  the  2  plane.  These  transformations  do  not  alter 
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of  that  expression.*  Substituting  (5)  into  (6): 

-Kir^=  -•  H-  f  ^ 

S  I  2rr 

This  is  not  the  imiqiie  solution  for  the  flow  in  this  problem. 

We  can  add  a  free  vortex  with  center  at  the  origin  and  arbitrary 
strength  *?*  .  We  eraress  its  velocity  potential  in  the  form 

—  ^  4- 

ZTT  C- 

whence  the  circle  still  remains  a  streamline  (with  ^/  =  O )  and 
the  velocity  at  infinity  is  not  changed.  This  additional  circu¬ 
lation  around  the  origin  will  enable  tis  to  take  account  of  a  pos¬ 
sible  Kutta  condition.  Thus  the  velocity  potential  (7)  can  be 
supplemented  and  reduced  as  follows: 

-  I, 


2w  O' 


To  evaluate  the  Blasius  integral  in  the  J  plane  we  note 

•  It  can  be  easily  shown  thet  the  imaginary  part  of  (6)  (the 
stream  function)  is  zero  on  the  circle,  thus  insuring  that  the 
circle  is  a  streamline.  Also  -W~{oo)  - 


And  therefore 

(13) 


'  [fs7  " 

-  Jtf)  -  TjT^l 

■  (^)  [  (I  j .)'  *  r  "j'l-iy " 

Ctry  f  e-^f  _  ^  c^ 

TT  L  5 


r  f'"’’  ^ 

e-cCP 

1 

"vllv^ 

\ 

1 

+  e*^c‘ 

jrj-f 

r  e-‘f 

b-i. 

?» 

.rC?-5.) 

r 

*?■(/- n 

-HLf-i—  _ I _ -1 

271'  If (j-^.;  f" 

5i 


+  mr_i —  j - ! _ 1 

2n‘Lfi-S.X/-i,')  Jfj-J.) 


’(r-5,Xf-|y)  f(f-U-) 


fNOT  REPRODUC 
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Equation  (10)  is  evaluated  by  deforming  the  contour  and 
evaluating  the  integral  around  the  singularities  of 

From  (13)  we  see  that  the  singularities  produced  by  are  all 

poles ,  located  at 


' 


and 

f. 


Since  the  vortices  are  located  at  ^  and  ,  these  points  are 
exterior  to  the  circle  and  residues  here  do  not  contribute 

to  the  integral.  On  the  other  hand|^|>C  guarantees  that  |  . |<  C.  , 

and  thus  the  residues  at  5^  =  (and  -sr-  )  must  be  accounted  for. 

The  singularities  of  (H))  all  lie  on  the  circle. 

They  are  poles  at  ^  s  t  C  and  branch  points  where  the  square  root  radical 
vanishes,  that  is  where 


or 


Hence  branch  points  are  located  where 


^  i  R 1 C  yc R' , 


that  is  at  the  four  points 
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where 


These  are  all  obviously  on  the  circle  since  in  each  case 
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SECTION  III.  Analytic  Evaluation  of  the  Contour  Integral 


If  a  complex  functionjp^f^,  has  a  pole  of  order  itj  at 
,  then  its  residue  there  is  (Ref.  15) 

(15)  (5  =  -^ 


Thus,  referring  to  Equ.  (10),  a  contour  around  a  pole  with  order 
>nof  the  Blusius  integral  located  at  ^  will  be 


(16)  ^  TTt 


C  f  (k 


where  ^  is  obtained  by  defining 


F(i)  ^^(i) 


in  Equ.  (15). 

The  poles  at  and  -C  are  simple  ones.  Thus  at  C, 


tne  integral  is 


T  n  t  P 


=  - 


^  ^  Cl-tC 


Similarly,  for  the  simple  pole  at  ^  =r  —  C- 


Thus  the  total  contribution  due  to  the  poles  at  C  and  -C  is 


2i 


CX>vio\;:sly,  terms  of  which  cire  either  real  or  odd  will  not 
contribute  to  this  expression.  Eliminating  these  terms  from 
Equ.  (13),  substituting  the  remainder  in  (18),  and  combining  terms 
gives,  after  considerable  algebraic  manipulation 


(19) 


I.  <  J|_  TrLi(c’-p 


ILL 


--  ^  I 


T  2c'Mrc‘- r/)‘ 

4-  f  -!-f ..  hi-.  ^  ,7  S 

4- 4-17  P,^  CP -  -I-  -ii 

_  zuf  -i— ,  4  XI } _ ! _ SI 

irU'-r  o'  Vc^-p/J 

+2xr_^4--LV _ L^s 

TT 

.  HJLl JL±LL^—-  -  — i 

IT  L  rc'-s/x 


«-  j 


s,  <  ^ 


>  <'’-3: 


+  XI  f-T^V  ^ 


..wi  kepkoducib 


The  contribution  to  the  total  moment  due  to  the  simple 
pole  at  JsO  is 

X  m/r[-  -  {T 


where  terms  of  (Equ.  (13))  have  been  included  which  contain  the 

factor  .  Obviously 

(20)  T./  =  ^ 

since  ^(c>)-o  .  (Equ.  (11).) 

To  find  the  contribution  of  the  second  order  poles  at  ^  ~ 
we  use  (16),  (15)  and  (17): 


(21)  To  A 


d.  r,. 


Ff-*"' 


C.  VP,(f f_ 

t 

1 

rr  \ 

s-s. 

•+  /: 

i  vp, 

f 

i 

j _ ! 

TT  \ 

/- 

r~  (Tr,-^lfP.  f  n, 


where  terms  of  IV  f$)  (Equ.  (13))  have  been  included  which  contain 
the  factor  -V  .  Since 

5 

and 

(see  Appendix  B),  we  have,  differentiating  (21),  letting  ^  s  0, 
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collecting  terms  and  eliminating  pure  imaginary  numbers: 


The  contribution  to  the  total  moment  due  to  the  third 


order  poles  at  <rc  is 

.  lT  P,  i  (P  '*  i  Pi  c  4  ^  y 

TT  /•  Jj 

v/here  terms  of  W^(%Jhave  been  included  which  contain  the  factor 
Hence 

«3,  i!  ^  A  (-?-/> 

since  l'*(o)~0  (Appendix  3). 

^  The  contribution  to  the  total  moment  due  to  the  fourth 
order  poles  at  ^  —  c  is 

r; /i  U...  y/<;^ 


where  terms  of 


•vY;) 


have  been  included  which  contain  the  factor 


Since  ^  fo)=  (Appendix  B)^ 

•r''  =  jli,  2irp^  cr><f  ■'“'’(r'’-  R^) 

J.O 

or 

(24)  r/  =  -  2f7; 


For  the  simple  poles  at  7  s  ,  the  contribution  to  the 
momentj again  using  (15),  (16'  and  (17)  is 


where  terms  containing  the  factor  - - in  the  expression  for 

have  been  inclxxied.  Th\is 

(25)  =  -A 

\  C' 


-  P.  W _ L 

ZTT  '  r  ^  -S 


where 


/(fi 


can  be  computed  using  Equ.  (11).  Similarly 


By  the  same  procedure 


so  that 


Fig.  9  Branch  points  and  branch  lines  of 
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The  remaining  singularities  where  the  contour  integral 
must  be  evaluated  are  at  the  branch  points  ^  and 

(see  Equ.  (14)).  We  arbitrarily  choose  as  branch  cuts  the  straight 
line  segments  joining  with  ^  (Fig.  9). 

The  function  which  contains  these  square  root  singu¬ 

larities  will  be  single  valued  and  continuous  as  long  as  w’e  move 
in  the  ^  plane  without  crossing  the  branch  lines.  Thus,  for  the 
integral  enclosing  points  '^^and  we  choose  the  contour  ABCDEFA, 
shown  in  Fig.  10,  where  CD  and  FA  actually  touch  the  branch  line 


Fig.  10 

but  have  been  shown  separated  for  visual  purposes.  £  will  be 
allowed  to  approach  zero.  Referring  to  equations  (10)  and  (11),  the 
integral  v;s  wish  to  evaluate  is 


\  £-  f  i  r  1 


or 


^gCDtf 


(28) 


To  evaluate  the  integrals  around  ABC  and  DEF  we  let  ^  ^  ^ 

and  t*'  respectively,  where  is  integrated  between 

0  and  27T’  in  the  first  case,  and  between  -  ft  and  JT in  the  second. 

If  the  integrand  is  denoted  by  then,  for  example, 

.-A-  ... 


It  nay  be  shown  in  general  that  if  /?  -  I  I 

^J(f-Cx')FU)|  =  o 

then  the  integral  vanishes  as  €^0,  Ihis  is  iisieed  the  case  for 
the  integrals  ABC  and  EEF  as  nay  be  verified  by  inspection.  (See 
Appendix  C  for  a  proof  of  the  theoren). 

Ifow,  considering  the  integrals  along  the  lilies  FA  and  DC, 
we  use  the  polar  coordinates 

(29b) 

where  '^^and  ^are  both  allowed  to  vary  between  o  and  irCwhich 
automatically  places  the  discontin’jity  along  the  designated  branch 
cut. )  On  FA 

=  O  and  ^ 

wltile  on  CD 

e^-2T  and  -n 


so  that 


V'jj-  varies  between  2R 


f  =  I 

The  only  change  in  the 
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Adding  (32)  to  (31),  and  after  algebraic  reduction,  we  obtain 


r  =  /  ^  / 


=  -  Hu  'f.pR' 


f  l')  =  0  ^  ‘  °' 


The  same  procedure  is  used  for  the  integral  around  the 
branch  cut  connecting  ^^with  and  produces 

(36)  I  ~""Kc.  J - rp - * - 


=  j- 


|r  ^ 


:37)  Ml')  =  • 
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It  should  be  noted  that  the  only  difference  between  the  inte¬ 
grands  of  (33)  and  (36)  is  in  the  sign  of  0-, 

We  may  combine  (33)  and  (36): 


(38) 

where 


(39) 


H  = 


If  we  note  that 


^  4' f  R*  (  1^  t- H  >p] 

^  o 

and  add  this  result  to  (38),  then  the  final  expression  is 

^  ^  ^ [ *■"( ’1)-'^  V >) )] 

J  .  •>  >£JJ^  ijg-  \  /C* 

with  given  by  (13),  H  (39),  ij  by  (34)  and  h  by  (35), 

Equation  (40)  is  a  definite  proper  integral  of  the  real  variable  K 
(but  with  a  complex  integrand)  and  must  have  a  finite  value. 

To  summarize,  the  rolling  moment  given  by  the  Blasius 
integral  is  the  svm  of  the  contour  integrals  around  the  singu¬ 
larities  inside  and  on  the  circle j^jsC  in  the  ^  plane.  That  is 

(41)  M.'  =  1/  t  x/  t  r;  ,  jj  ,  X 


+  1;./-  f  I,,,, 


c 


Vi. 


•t-  J4.  ,  )  T 


NOT  REPRODUCIBI  ^ 
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The  moment  due  to  fin  cant  must  be  added  to  this  result.*  The 
equation  for  each  of  the  integrals  in  (41)  is  listed  in  the  table 
below: 


Integral 

Equ.  nianbar 

I.'  ’ 

(19) 

x;  2  c 

(20) 

t: 

(22) 

o 

III 

(23) 

T 

(24) 

T  ' 

^  <■'/%. 

(25) 

Vi. 

(26) 

Jo 7^  ±  Xcyf^ 

(27) 

r.  4  X. 

1 

(40) 

The  rolling  moment,  Equ.  (41),  can  best  be  evaluated  for 
various  values  of  the  parameters  by  use  of  a  large  scale  digital 
computer.  However,  we  can  obtain  iresults  directly  for  special 
limiting  cases.  For  example,  consider  the  two  dimensional  flow, 
in  the  absence  of  external  vortices,  past  a  flat  plate  inclined 
at  a  roll  angle  (angle  of  attack)  ^  with  respect  to  the  free 


•See  Section  II  and  reference  3. 
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stream  velcx:ity,  TJ  .  Obviously,  since  Fig  (4b),  Equ.  (40) 


indicates  that 


^  Ir  <-  =  °- 

Jr.  ?Ti  for,  Jn:  n  n 

Similarly,  no  external  vortices  in  the  flow  means  that  «  i  —lx 

consequently,  all  the  ether  integrals  in  Equ.  (41)  are  clearly 

. —  H 

zero,  with  the  exception  of  J_^  .  Thus,  from  (24) 

-1/  =  -  IP  TT  Cr’-a.^^'v-  S.  f  • 


But  for  R  =  C 


Ql  ^  C  = 


Hence  for  a  flat  plate  with  no  vortices 

/I  ^  TT'tr^  4  ^  ^  cp 

(42)  /  ^ 

=  -  ^  t?  ^  ^  ^ 

which  is  the  well  known  moment  about  the  center  chord  in  ideal 
flow.  In  this  case  the  lift  L  acts  at  the  quarter  chord  (aero¬ 
dynamic  center)  so  that 

l'4  =- 


whence 


L  =  _  ±kL 


or  using  (42) 


which  becomes,  for  small  angles 

C.,  = 
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r  I 


I 


the  classical  result  (see  Ref.  16). 

The  other  limiting  case  we  will  consider  is  that  of  a 
f inless  cylinder.  Here  ^  =  R,  (Fig.  4)  so  that 

C  =  ^  snd  Cl  =  O , 

If  there  are  no  external  vortices  in  the  flew,  =  O  and 


=  I/  +  r 


•h  T 


Since  CL  =  0,  Equ.  (24)  shoivs  that 


T  =  0 


Also,  from  (34) 


and  from  (35) 


7|  (rl  =■>] 

=  CV  z'yR-r). 

Substituting  these  values  into  (39)  gives 


H  = 


-  (R-f) 

1 

[(R-k)'- 

-R] 

J 

( 

\.-y‘ (2R-i)‘+ 

I'-fR" 

■] 

which  is  pure  imaginary. 
Moreover 


yjj  -  W^(  >j) 


=  Cv'  ( ’j) 

r.  -  1  QyyL  ^*^'>1). 
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Thus,  Equ.  (40)  indi.cates  that  the  integrand  is  the  product  of  two 
imaginary  variables,  and  is  therefore  real.  Hence 

=  /it  4-  ‘'P  ^  /  (Real  variable) 

4 


ir,  So:  f, 


=  fh  i  (Real  number) 


=  O 


and  so  the  rolling  moment,  Mj  ,  is  zero  for  a  finless  cylinder, 
the  expected  result  in  ’.deal  (frictionless)  flow.* 


•It  is  also  expected  that  the  rolling  moment  would  be  zero  even  in 
the  presence  of  vortices,  i.e.,  P,  and  not  zero.  In  this  case, 
since  C*i^,  ClsO  and  fay-  it  can  be  shown  that  Equ.  (19)  and  (27) 
both  vanish  while  Equ.  (22)  cancels  (25)  and  (26)  so  that,  indeed. 
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SECTICTJ  '^V  Conclusion  and  Plans  foe  Future  Work 


An  analytical  expression  for  the  sectional  moment  rolling 
coefficient,  A/J,  has  been  found  as  a  function  of  the  following 
variables : 

There  are,  in  fact,  twelve  independent  Vaj-iabies  since  ^  and  ^  ^ 
are  complex  numbers,  each  consioting  of  two  components.  We  may 
eliminate  three  parameters  if  we  express  the  variables  in  dimen¬ 
sionless  form: 


(43) 


=  C^.f  -L  -L  Cp  Or  q> -2—- 


where  "S 


P/zttR  cr 


is  defined  as  the  vortex  strength  coefficient.  It  should  be  noted 
that 

V 

where  the  extremes  B.  s  is  for  a  flat  plate  and  *4'  =  |  is  for  a 
f inless  cylinder.  The  numerical  data  obtained  from  a  computer  will 
be  most  conveniently  plotted  using  (43).* 

As  an  example,  we  suppose  that  r  is  determined  by  a  Kutta 
condition  and  so  need  not  be  specified  independently.  For  vortices 
fixed  at  and  ,  we  may  plot  series  of  graphs  with  as  ordinate , 
Cf  as  abscissa  and  ^  as  parameter,  where  each  graph  is  for  a  pair 
of  values  Cj,  and  Cp^.  Thus  for  the  case  c-  we  have,  from 


*  The  variables  R  and  (r  are  defined  in  the  physical  plane  (see  Fig. 
4).  If  the  location  of  the  vortices  are  at  "2, and  ?jin  the  physical 
plane,  their  locations  in  the  ^  plane  may  be  determined  from  Equ. 
(4a).  It  will  also  be  recalled  t'natV^cp  and  P  suffer  no  change 

XU  uKw  \,.,kOii5A.  wxitiacxv/il* 


(42),  for  a  flat  plate 

(f=°) 

while  for  a  finless  cylinder,  there  is  no  rolling  nesnent  in  friction¬ 
less  flew  so 


Fig.  11 

The  analysis  in  this  report  assumes  a  two  dimensional  fin. 
However  if  the  fin  is  not  rectangular,  but  has  varying  sp)an  (r  ,  we 
may  integrate  the  sectional  rolling  moment  H*'  in  steady  flo-/  to 
find  the  total  moment 


where  fid  is  a  function  of  iu)  and  L.  Is  the  fin  length.  (See  Fig 
12).  Then  in  Equ.  (43)  we  may  replace  the  dimensionless  variable  (L|i^«by 
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R 

L 


Fig.  12 


In  conclusion  it  is  to  be  noted  that  this  report  is  the 
first  in  a  continuing  research  effort.  The  immediate  future  task 
is  to  obtain  namerical  results  for  the  present  analytical  study. 
As  another  important  undertaking  it  will  be  useful  to  extend  this 
work  to  missiles  with  a  cruciform  tail  body  configuration.  Also 
compressible  and  non-steady  flow  effects  should  be  considered. 
Finally  meaningful  experiments  should  be  devised  and  carried  out 
to  supplement  the  analysis. 
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APPENDIX  A 


It  will  be  proved  here  that  a  vortex  undergoes  no  change 
in  strength,  P  ,  under  a  conformal  transformation  from  the  2  to 
planes.  By  definition: 

^  f  u  J'X  +  V-  ) 

where  the  integral  is  evali;ated  comtercloclcwise  in  a  cxirve  en¬ 
closing  the  vortex. •  Now,  in  the  complex  2  plane 


I*  (  ^ }<?!  2. 


-  (  U  —C  v)(c^y  4- 

(  C(  dy  i  V  j 


so  that 


plane 


Similarly,  in  the  ^  plane 


But  by  Equ.  (9) 


w(  t)  d  2,  —  f  i )  d  ^  ^ 


Hence  it  follows  that 


i  plane 


•  The  negative  sign  is  due  to  the  fact  that  a  positive  circulation 
is,  by  convention,  defined  to  be  in  the  clockwise  direction. 
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APPENDIX  B 

For  reference  pvirposes,  the  derivatives  of  /'ff)  are  listed 
here.  From  Equ.  (11) 


/«’ "  *  - ; 


where 


Hence 


i  r‘  - 


ve'j  '  +  Cf'-'ic-'; 


where 


Also 


/  •a;  .  ph  j0^  ;  n; .  jfe.  7  >!/ 


and 


/7  /« 


f 


-  3 


elii... 


It  is  easy  to  show  that 


(  = 


=  ^ 


•(O)  =  O 


'  NOT  REPRODUC 
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APPEl^IX  C 

Suppose  we  wish  to  integrate  F(^^in  a  circle  of  increasingly 
infinitesimal  radius  around  a  branch  point,  such  as  in  Pig.  10. 

It  will  be  proved  that  this  integral  vanishes  if 


Hence 


/n 

\\^e-^F 

^->0  J c 

1  i-n 

o 

l(^^r) 

r  ( 5}|  ci-o- 

1(5 

i.)  F 

'  ITT 

Tx 

M4X 

where  y  is  that  point  on  the  circular  contour  which  gives  the 
maximum  value  indicated.  Now  if 

^  1  =  o 
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for  any  mode  of  approach  of  ^  to  ,  then  it  is  also  true 

along  the  path  of  maximum  absolute  values.  Hence 

^  pfi)  e/'?  =  o 

and  the  theorem  is  proved. 
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